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subsequent  decay  of  the 


seasonal  tl'.ernocline  is  acconpanied  by  iiany  rand  on  fluctuations.  The 
physical  processes  causing  the  decay  and  the  fluctuations  in  the  nixed- 
layer  depth  are  cxai.uned  and  qualitatively  evaluated, 

Tne  Tukey  spectrun  analysis  progranned  for  the  CDG  1604  electronic 
digital  computer  was  used  in  analysis  of  oscillations  at  the  bottom  of 
the  mixed  layer  to  determine  the  distribution  of  \;ave  energy  for  series 
of  data  taken  at  hourly  intervals.  -An  energy  peak  centered  neai-  12  hours 
was  obsei’\’’od  to  prcdoi-ainate,  and  this  energy  is  equiva.lent  to  that  of  a 
12-hour  internal  \;ave  \;ith  height  of  27,0  foot.  The  computer  \;as  also 
used  to  select  the  significant  meteorological  and  oceanographic  parameters 
which  could  be  used  to  predict  a  change  in  the  mixed “layer  depth  by  iise 
of  the  BUD  07  multiple  regression  proqprrm.  The  most  significant  paa’ameter 
was  sea  surface  temperature.  The  best  correlation  coefficient  for  this 
parameter  occurred  for  lags  of  aero  to  12  hours. 

Uie  results  of  the  analysi  s  support  the  theory  that  convection  is 
the  physical  process  which  causes  the  seasonal  decay  of  the  thermo- 
cline  dtiring  the  cooling  season  and  that  short-term  fluctuations  of  the 
mixed-layer  depth  are  due  primarily  to  internal  \;avec. 
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1 . 

The  vertical  tenporature  structure  of  the  near- surf ace  layer  in  the 
open  ocean  is  a  conplicated  phenomenon  v;hich  varies  in  space  and  tine. 

Certain  aspects  of  the  temperature  structure  arc  of  great  significance 
in  determining  sonar  effectiveness.  Thus,  both  the  military  and  the  com¬ 
mercial  fishing  industries  are  interested  in  prediction  of  the  thermal 
structure. 

Various  processes  continually  modify  the  thermal  structure  to  produce 
an  essentially  homogeneous  (mixed)  surface  layer,  underlain  by  a  thin 
layer  of  rapidly  increasing  density,  v;hich  varies  in  thickness  through- 
o'af  the  year,  Tlic  aixed«»layer  depth  (referred  to  as  the  I-ILD)  is  the  depth 
belov/  the  v;atcr  surface  to  which  mixing  has  established  isothermal  condi¬ 
tions,  The  lov/er  boundary  of  the  IILD  is  the  top  of  the  thormoclino,  a 
thin  layer  of  largo,  negative  vortical  temperature  gradient,  usually  associa¬ 
ted  with  the  layer  of  increasing  density,  A  complex  interaction  of  meteoro¬ 
logical,  oceanographic,  and  (indirectly)  astronomical  factors  provide  the 
driving  forces  for  changes  in  the  ocean  thermal  structure. 

Random  and  periodic  fluctuations  of  the  IILD  have  been  examined  by 
many  investigators.  It  is  knovm,  as  a  result,  that  the  iILD  may  vary  from 
a  fcv7  inches  in  the  heating  season  to  several  hundred  feet  in  the  cooling 
season.  At  times  it  may  fluctuate  tens  of  feet  in  a  matter  of  minutes. 

The  objective  of  this  research  is  to  identify  more  clearly  significant 
processes  which  affect  the  IILD  during  the  cooling  season  and  to  relate  the 
IILD  to  meteorological  and  oceanographic  parameters.  The  processes  capable 
of  affecting  the  vertical  temperature  distribution  are  discussed. in  the 
appendix.  Certain  of  these  processes  will  be  considered  in  this  paper: 

(l)  thermohaline  convection,  (2)  dynamic  convection  due  to  wind  and  vavo 
action,  and  (3)  unidentified  processes  which  generate  internal  waves. 
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Up  to  now,  unly  yaalit  t-VJ  uoocri^jtiono  of  the  ..axod-loycr  behavior 
durin{:;  the  cooling  season  are  foiind  in  the  litcratui-e.  These  usually 
state  that  increased  wind  stirring  and  surface  cooling  act  together  during 
the  cooling  season  to  drive  the  seasonal  thcrraoclinc  do\.'n  to  its  riaxinun 
annual  depth.  Our  findings  will  shovr  that,  at  Station  Paoa,  thernohaline 
convection  is  the  nest  important  factor  and  that  v;ind  nixing  (dynamic  con¬ 
vection)  is  actually  insignificant  during  the  cooling  season. 

Met e Orel ogical  and  oceanographic  pararacters  contributing  significantly 
to  the  change  in  the  MLD  during  the  cooling  season  arc  deterninedj  and  from 
those  the  physical  processes  of  importance  arc  inferred.  The  investiga¬ 
tion  was  carried  out  with  the  aid  of  statistical  correlation  and  wavo- 
analysis  programs,  the  BIID  07  [l^  and  the  lUkcy  spectrum  analysis  [153, 
utilizing  the  GDC  I6O4  digital  computer.  Regression  equations  and  correla¬ 
tion  coefficients  arc  obtained  by  the  BIID  07  program.  One  part  of  the 
Tukey  spectrum  analysis  was  used  to  determine  cross  spectra  of  the  IDD 
fluctuations  and  the  sea-surface  temperatures.  This  program  was  also  used 

t 

in  the  analysis  of  internal  v;avcs. 

The  investigation  is  ba.scd  on  oceanographic  and  weather  da.ta  obtained 
by  the  Fisheries  Research  Board  of  Canada,  Pacific  Occanogi’aphic  Group, 
for  the  months  of  October,  flovcmber,  and  December  during  the  years  1956, 
1957  and  1958  [6.  7]. 
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2 .  '3ackf,Tound 

The  attempt  to  identify  tlic  factors  v/hich  bring  about  the  ho\irly 
fliietuations  in  the  lILD  v;as  based  on  the  authors'  earlier  endeavor  to 
derive  an  empirical  method  to  predict  the  average  MI.D  on  a  daily  basis 
during  the  cooling  season.  It  became  apparent  then  that  the  hourly  fluc¬ 
tuations  of  the  I  ID  could  be  of  equal  or  larger  magnitude  than  the  average 
deepening  of  tlio  tlLn  during  one  or  more  days;  thus,  to  use  one  or  an 
average  of  tvra  or  more  observed  I'XD  values  per  day  as  representative  of 
the  depth  for  any  particular  hour  of  that  day  is  misleading  Q  J  .  The 
probable  cause  of  these  fluctuations  is  internal  v;avcs. 

Students  at  the  U.  S.  Ilaval  Postgraduate  School  during  1961  and  1962 
derived  empirical  methods  to  predict  the  M'D  at  Ocean  Station  Papa  dui’ing 
each  of  the  12  months  of  the  year.  In  addition,  the  onset  of  the  seasonal 
thcrmoclinc  at  Ocean  Station  Papa,  uhicli  talces  place  in  April  or  May,  v:as 
studied  by  Clark  in  1961  [2]*  Clark  c.uccessfull.y  related  the  thermocline 
onset  to  upper-air  parameters.  During  the  heating  season,  v;hen  the  sea¬ 
sonal  thormoclino  gradually  increases  in  magnitude,  vjind  mixing  processes 
have  been  clearly  established  as  the  dominant  factor  in  dotoiraining  the 
I  JiD,  by  Tabata  [22^  end  Geary  [lo]]  in  1961.  The  cooling  season  remains  to 
be  studied.  V/e  have  made  this  attempt. 


3.  Dcccrijition  of  Ocoan  St-ition  Papa 

The  geographic  location  of  Ocean  Station  Papa  is  50N,  145W  (figairo  l). 
It  is  approximately  400  miles  north  of  the  boundary  botveon  the  eastern 
subarctic  Xv'ater  roass  and  the  subtropic  uator  rnss  [8]  .  The  vrater  depth 

at  Papa's  location  is  2000  fathoms;  therefore,  the  thermal  structure  is 
in  no  v/ay  influenced  by  bottom  topography.  Also,  its  location  is  approxi¬ 
mately  in  the  center  of  the  A].askan  Current  G’/ral  Q  9]  (figure  1). 
Because  of  this  location,  the  effects  of  advcction  arc  minimized  apd  will 
not  be  considered. 

The  ueather  ship  normally  maintains  an  on-station  position  \;ithin 
a  ter— mile  square  centered  on  50H,  145W,  In  order  to  avoid  rough  weather, 
the  ship  may  steam  anyv/hore  \;ithin  the  210-milc  square  grid  (figure  2), 
Sach  bathythermograph  (3T)  sounding  is  identified  by  a  tvro-lotter  code. 

If  the  ship  is  on-station,  the  letter  group  OS  indicates  the  position. 
The  ship  position  for  the  greater  part  of  the  period  of  our  concern  v;as 
v/ithin  a  30-r,ilc  square  centered  at  OS.  .Iny  data  outside  of  this  30-milo 
square  were  not  considered  in  this  study. 
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4. 


oi  To;..' •'r  •.  1..-’.  S'llinity  oiructui-c  dvi’inj  tl'.o  Cooling  Sc.inon 
x'^iguvo  3  ro,,rooontr;  tho  docny  of  the  thorr-ioclino  during  the  yoor 
1956.  In  A’agunt,  the  rurfacc  Inyoro  have  boon  heated  to  a  naxinura,  approxi¬ 
mately  37F.  Due  prinai'ily  to  wind  mixing,  there  is  a  sliallou  nixed  layer 
approxirately  ICO  foot  in  depth  with  a  sharp  underlying  thernoclinc.  With¬ 
in  the  thcrnocline,  the  water  temperature  drops  to  4OF  at  a  depth  of  200 
feet.  Progressing  into  the  fall  and  vdnter,  the  M-D  continually  increases 
and  the  magnitude  of  the  thernocline  docrcasoc.  By  the  end  of  December, 
tlae  MTjD  is  nearly  330  feet  in  dcptl:^  and  the  sea-surface  tonpcratui’o  has 
uccroasod  to  approximately  42F.  Due  to  the  presence  of  an  intense  halo- 
cline,  the  MIjD  I'cach.os  a  limiting  depth  of  330  feet  in  Jantiarj^  and  iso¬ 
thermal  conditions'  prevail  down  to  this  depth  for  the  remainder  of  the 
winter. 

The  basic  salinity  striicturc  for  the  eastern  subarctic  Pacific  v/atcr 
mass  is  shown  in  figure  4.  Tliroc  zones  arc  indicated.  The  upper  zone 
extends  from  the  sea  surface  to  a  depth  of  320  foot,  and  is  characterized 
by  rolatively  lo\;-salinity  v;ator  (32, '59^  f2l]  la  the  fall,  the  \;ppcr 
zone  is  isoh.aline  and  is  narlceu  by  the  presence  of  the  seasonal  thermo- 
clinc.  Dccaixso  of  the  Isohallnc  conditions  in  the  upper  330  feet,  vre 
assume  that  significant  variation  in  density  is  governed  by  changes  in 
temperature  only. 

The  haloclinc  represonts  a  transition  zono  botv/oon  the  upper  and 
lower  zones.  Kero,  the  salinity  increases  markedly  with  depth  (as  much 
as  one  part  per  thousand  v.’ithin  .on  interval  of  330  feet).  The  haloclinc 


represents  a  limiting  depth  for  the  seasonal  thormoclinc.  It  is  similar 
to  having  an  ocean  bottom  at  this  depth  as  far  as  convective  nixing  is 
concerned.  Tlxe  density  increase  is  so  great  through  this  layer  that  mic'ing 

cannot  take  place  through  it  ^20] .  Within  the  lower  zone  0^50  feet) 
the  salin  i, -"adually  incroac 


cs  xdth  depth  to  the  ocean  bottom, 
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*  C  .j  j  oc  t 

Tho  purpose  of  this  statistical  study  of  the  MLD  is  three-fold. 

All  attnnpt  is  mde:  (l)  to  dctorirdne  the  relative  inportanec  of  thcrr.o- 
halinc  conveetive  mixing  and  wind-inducod  mixing  on  tho  deepening  of  tho 
thormoclino  during  the  cooling  season;  (2/  to  determine  tho  lag  betv;oen 
a  signifioant  meteorological  or  oceanographic  event  and  a  subsequent  cliango 
in  the  MLD;  and  (3)  to  examine  tho  inportanoo  of  the  seemingly  random 
short-term  fluctuations  of  the  MTjD  and  to  provide  a  description  of  those 
internal  v;avcs  present. 

First,  the  parameter  or  parameters  are  identified  which  have  the 
largest  linear  correlation  viith  the  deepening  of  the  MiD.  Next,  a  further 
investigation  has  been  undertalcen  of  these  significant  parameters  to  deter¬ 
mine  the  lag  in  the  mixing  process,  that  is,  the  time  between  the  onset 
of  a  change  in  the  parameter  and  a  subsequent  deepening  of  the  MLD. 

Finally,  the  observed  variations  in  the  M.D  v/hich  eannot  bo  aeoounted  for 
by  either  mixing  process  arc  investigated, 

5.b,  Method  of  Investigation 

The  parameters  to  bo  investigated  are  sub-divided  into  three  groups. 
Those  related  to  convective  stirring  arc:  (l)  temperature  of  the  sea, 
at  tho  time  of  the  sounding  and  v;ith  a  2-hour  lag;  (2)  air  temperature, 
at  the  time  of  tho  sounding  and  v;ith  a  2-hour  lag;  (3)  devupoint  tempera¬ 
ture  and  v;ct-bulb  temperature,  at  the  time  of  the  sounding;  (4)  relative 
humidity,  at  2-hour  and  5-hour  lags;  and  (5)  low-eloud  cover,  with  a  2- 
hour  lag. 

Parameters  chosen  to  represent  v;lnd  and  v/ave  mixing  are:  (1)  v;ind 
speed,  at  2-hour  and  5-hour  lags;  and  (2)  height  of  significant  waves, 
at  2-hour  and  5-hour  lags. 


Parameters  chosen  ’.,’hich  verc  chov;n  to  have  a  correlation  v;ith  the 
th-crmal  structure  by  others  in  previous  studies,  but  did  not  conveniently 
fall  into  the  abo'^’c  categories,  are;  (l)  500-mb  height,  in  goopotential 
meters;  (2)  500-rab  temperature  and  500-nb  relative  hixmidity,  at  a  2-hour 
lag;  and  (3)  surface-pressiu’o  change,  v;ith  a  2-hour  and  5-hour  lag. 

Those  data  were  obtained  from  synoptic  3-hour ly  reports,  tv,’ice-daily 
constant— pres sure  records,  and  the  BT  traces.  The  available  BT‘s  for  the 
period  of  investigation  wore  tabulated  in  tvjo  groups,  0200S  and  17003, 
for  each  year  to  eliminate  any  possible  yearly  anomalies  or  diurnal  trends. 

An  investigation  by  Soripps  Institution  of  Oceanography  of  the  annual  MLD 
fluctuation  in  the  northeastern  Pacific  Ocean  shows  that  the  MLD  most  jfre- 
quontly  observed  in  a  given  month  often  varies  oonsidorably  from  year  to 
year  [38].  The  available  literature  on  diurnal  trends  during  the  cooling 
season  is  not  so  clear.  A  preliminary  investigation  of  the  data  showed 
no  marked  diiornal  trend.  This  indication  v/as  supported  by  Dr.  Tully  of 
the  Pacific  Oceanographic  Group,  in  a  conversation  with  the  authors,  in 
v/hich  he  stated  thau  during  the  cooling  season  the  MLD  is  too  deep  to  have 
a  clearly  defined  diurnal  fluctuation, 

With  the  aid  of  a  multiple-regression  and  correlation  analysis,  the 
BIMD  07  programmed  for  the  GDC  1604  high-speed  eomputer,  the  linear- 
correlation  coefficients  were  obtained  for  the  prc-solcctcd  group  of  me¬ 
teorological  and  oceanographic  parameters  with  the  MLD.  The  BIMD  07  pro¬ 
gram  was  used  because  it  can  select  different  sub-samples  of  data  that 
are  obtained  from  the  sane  population.  Thus,  the  program  would  perform 
linear-coirolation  analysis  and  linear  mltiple  regression  on  the  data 
for  02003  and  17003,  separately  and  in  combination.  In  this  v.'ay,  a  close 
examnation  of  possible  diurnal  trends  could  be  made  for  each  of  the  throe 
years  studied,  1956  through  1958.  The  analysis  provided  regression  equations 
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for  nine  difforent  data  groups  for  the  periods  shown  in  Table  I. 

The  BII'D  07  program  yields  best  results  v;hcn  the  nui.iber  of  indepen¬ 
dent  variables  is  not  greater  than  one-half  the  number  of  values  of  the 
dependent  variable.  Since  the  smallest  number  of  values  for  the  depen¬ 
dent  variable  is  36  (at  02003,  1958) only  18  independent  variables  were 
investigated.  A  set  of  18  secondary  parameters  \;as  tabulated  in  the  event 
that  the  analysis  of  the  18  primary  parameters  proved  inconclusive. 

Tlicrc  was  no  attempt  by  the  authors  to  eliminate  possible  duplication 
of  parameters  that  arc  indicative  of  the  same  process.  For  example,  dc\-^- 
point  temperature,  wet-bulb  temperattirc,  relative  humidity,  and  low-cloud 
cover  are  all  related  to  evaporationj  but  it  was  left  to  the  BIl-ID  07  program 
to  discern  which -ones,  if  any,  were  significant  in  contributing  to  the 
prediction  of  the  deepening  of  the  liLD. 

The  BIM)  07  is  unable  to  produce  a  multiple  linear-regression  equation 
if  any  independent  variable  is  a  linear  combination  of  tv/o  or  more  inde¬ 
pendent  variables.  Therefore,  temperature  of  the  air  minus  temperature 
of  the  sea,  \jhich  is  often  used  as  an  indicator  of  thcrmohalinc  convection, 
could  not  be  investigated  in  conjunction  vzith  the  air  temperature  and  sea 
temperature.  Consequently,  three  additional  problems  were  computed  using 
the  air-minus-scc  fcenperature  parameter  in  place  of  sea-temperature  para¬ 
meters.  See  Table  II  for  the  tabulated  results  of  the  12  problems. 

Tn  analysis  of  all  the  data,  certain  facts  and  tendencies  ore  evident; 
these  arc  discussed  in  the  following  two  paragraphs. 

The  only  parameters  that  shov;  a  consistently  high  linear  correlation 
with  the  MI.D  in  every  sample  or  sub-sample  arc  the  sea-surface  tempera¬ 
tures,  at  the  time  of  g  sounding  and  v/ith  a  2-hour  lag.  They  also  accounted 

"*(Tliis  small  num.ber  results  because  BT  data  for  October,  195§  arc  miss¬ 
ing;  only  a  small  nmber  of  observations  was  available  in  1957  as  vjcll,  since 
December,  195*^  observations  were  not  taken. ) 
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for  alnoGo  all  of  tiio  varianco  explained  by  the  regreGnlon  equation  using 
all  18  paraiaetcrs.  Figure  4  shows  the  close  relationship  bctv/cen  the 
sea- surface  tenporature  and  the  MLD  from  October  through  December,  1958, 
Likcvfisc,  Zcttcl  |253  suggested  that  sea- surface  tomnerature  might  be  a 
good  "indicator"  for  the  IXD. 

V/ind  speed  and  height  of  the  significant  waves,  the  parameters  chosen 
as  indicators  of  wind  stirring,  have  either  very  lovr  negative  or  lov;  posi¬ 
tive  correlation  coeffieients. 

In  order  to  determine  the  time  required  for  an  increase  in  convoetivo 
activity  due  to  a  reduction  of  sea-surface  temperature  to  effect  a  sub¬ 
sequent  deepening  of  the  MjD,  the  llLD’s  for  17002,  tabulated 

v;ith  the  sea  temperature  at  lags  of  zero,  two,  five,  eight,  eleven,  four¬ 
teen,  seventeen,  tv/enty-onc,  and  tv;enty-three  hours  as  the  independent 
variables.  A  simple  linear-eorrolation  prograra  was  used  to  determine 
that  the  lag  yielding  the  best  correlation  was  two  hours  (Table  III). 

Also,  two  one-hourly  series  of  IILD  for  79  hours  in  October,  1957,  and 
for  51  hours  in  November,  1958,  v;cre  tabulated  v;ith  their  eorresponding 
sea-surface  temperatures.  The  energy-density  spectrum  and  the  eross-spcc- 
ti*um  of  both  the  IITjD  and  the  sea- surface -temperature  time  series  i:as  ob¬ 
tained  from  the  Tukey  analysis.  Maximum  lags  of  39  hours  and  24  hours 
v;ere  used  for  the  79-hour  and  51-hour  series  respectively.  The  fv/o  co- 
spectio  for  the  79-hour  and  the  51-hour  series  show  largo  values  for 
h’s  around  8  hours  and  4  hours  respectively  (figure  5)» 

From  the  results  of  the  simple  linear  correlation  and  the  Tukey  analysis,^ 
it  is  soon  that  all  lags  from  two  to  six  hours  show  high  correlation.  8c- 
yond  12  hours  the  correlation  rapidly  decroaser, 

A  linear-regression  equation  was  derived  for  17002,  1956,  using  sea 
temperature  v;ith  a  2-hour  lag  as  the  single  independent  variable  (Table  IV  ). 


1  on  cross-correlation,  not  shown 
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Tho  fact  that  the  results  of  this  equation  as  measured  by  the  standard 
error  conuarc  favorably  with  tho  multiple- regression  equation  derived  from 
13  independent  variables  for  the  same  year  supports  the  initial  analysis 
of  the  problem.  If  a  regression  equation  vjsrc  to  be  used  to  forecast  1I.D, 
one  using  only  the  sea- surface  temperature  taken  with  a  lag  of  up  to  12 
hours  would  be  much  more  reliable  than  an  equation  with  more  parameters. 
This  is  true  simply  because  it  would  have  a  larger  number  of  degrees  of 
freedom.  Note  that  all  of  the  regression  equations  have  a  standard  error 
of  from  15  to  30  feet,  with  an  average  for  all  problems  of  22.15  feet 
(Table  V,  Table  VI). 

There  are  several  possible  explanations  for  this  standard  error  of 
more  than  20  feet:  (l)  ft  is  possible  that  not  enough  parameters  v/crc 
used,  and  that  by  substituting  for  primary  parameters  v/hich  cho\/ed  a  very 
poor  correlation,  some  of  tho  secondary  parameters  or  primary  parameters 
with  different  lags,  the  standard  error  of  the  new  regression  equation 
would  be  reduced;  or  (2)  there  is  a  strong  possibility  that  a  non-linear 
relationship  between  a  parameter  and  the  MLD  could  account  for  a  large 
part  of  the  standard  error  (hov/cver,  the  BIMD  07  program  cannot  detect 
this);  or  (3)  the  error  may  be  due  to  internal  waves.  Of  these  possi¬ 
bilities,  (3)  seems  tho  most  likely  reason  to  explain  tho  observed  stan¬ 
dard  error. 
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6.  ilnalysis  of  Intornal  VJavoc 

The  base  of  the  jaLxed  layer  in  the  ocean  is  cubjcct  to  nany  vertical 
oscillations.  Figure  7  sho\Js  a  series  of  daily  MLD's  for  the  period  Octo¬ 
ber  through  Dcccnbor,  1956.  Each  plotted  point  on  the  curve  is  an  average 
daily  value  of  the  IXD.  Tv;o  individual  soundings  arc  normally  used;  how¬ 
ever,  up  to  24  soundings  wore  available  on  several  days.  It  is  observed 
that  the  MIjD  undergoes  a  seasonal  downviard  trend.  During  the  October 
through  December  period,  the  seasonal  trend  of  the  MjD  finds  it  increas¬ 
ing  from  approxiriately  130  feet  to  330  feet.  Tliis  is  an  average  deepen¬ 
ing  of  2.17  feet  per  day.  Also,  largo  vertical  fluctuations  are  super¬ 
imposed  on  this  dovmuard  trend,  v/hich  become  larg«‘"  in  magnitude  as  the 
nixed -layer  thickness  increases.  Note  in  figure  7  that  the  vertical  fluc- 
tixations  of  the  ILD  can  be  ten  times  greater  than  the  average  daily  deepen¬ 
ing.  In  addition,  it  is  seen  (figure  8)  that  hourly  fluctuations  can  bo 
of  equal  or  greater  magnitude  than  the  average  daily  deepening  during  the 
fall. 

Figure  8  represents  a  one-hourly  tine  series  for  51  consecutive  hours 
during  November,  1958  (series  l).  A  second  one-hourly  time  series  rep¬ 
resents  79  consecutive  observations  taken  during  October,  1957  (series  2). 
The  magnitude  of  the  IILD  fluctuations  is  about  the  same,  whether  it  ir 
computed  from  hourly,  twice-daily,  or  daily  observations.  The  average 
fluctuation  per  hour  of  the  MLD  for  series  1  is  18  feet.  The  average 
fluctuation  between  the  consecutive  02003  and  17003  soundings  for  the  month 
of  November,  1958^ is  21  feet.  Also,  the  averaged  daily  fluctuations  during 
November  is  25  feet.  The  nurabers  for  series  2  arc  similar.  A  December, 
1956, one-hourly  tine  series  (figure  9)  Was  available  but  was  too  short 
for  detailed  analysis.  Hov;ever,  it  illustrates  the  extreme  fluctuation 
of  the  I'iLD  v;hich  talces  place  over  a  short  time  interval.  Over  a  3-hour 
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,)ci*iod,  tl. .  1 LD  incr'’'ircd  fror.i  270  To  t  to  370  foot,  and  then  decreased 
attain  to  235  foot. 

Since  nixing  processes  can  account  for  only  an  increase  in  the  IXD, 
and  since  observed  fluctuations  are  in  both  directions  (up  and  do\;n), 
and  arc  a  consistent  feature  of  the  data,  these  observed  fluctuations 
arc  considered  to  bo  due  to  internal  vaves. 

Internal  v/aves  nay  occur  \jithin  stratified  water  and  in  water  in 
\;hich  the  density  increases  \iith  depth.  The  largest  vertical  displace- 
nents  of  tlie  iratcr  particles  arc  to  be  found  v/ithin  the  boundary  between 
layers  of  different  density.  The  displacement  ar:iplitudo  diminishes  above 
and  bclov/  this  boundary,  approaching  zero  at  the  sea  surface  and  the  ocean 
bottom.  The  density  difference  betv/een  the  waters  above  and  below  the 
seasonal  thcrmocline  is,  of  coiuse,  much  less  than  the  density  difference 
between  the  air  and  the  v;ater  at  the  sea  surface.  Therefore,  a  boundary 
v/ithin  the  ocean  can  be  much  more  easily  displaced  than  the  surface  of 
the  sea  CL7]  Ufford  [24]  states  that  surface  v/aves  require  30,000  times 
more  energy  to  start  and  maintain  than  is  required  for  internal  v/aves  of 
the  same  amplitude.  Energy  for  internal  v/aves  could  come  from  surface 
disturbances  or  from  currents  v/ithin  the  ocean. 

Besides  rapid  vertical  oscillations  in  the  IIjD,  it  v;as  observed  that 
the  thickness  (or  vertical  extent)  of  the  thcrmocline  layer  bclov/  the 
lILD  may  vary  markedly  v/ith  time.  Figure  10  shov/s  the  thcrmocline  thick¬ 
ness  for  scries  1 .  Variations  in  the  thcrmocline  thickness  indicate  that 
intcrnal-v/ave  characteristics  of  phase  and  amplitude  vary  v/ith  depth, 

Mso,  the  temperature  gradient  in  the  thcrmocline  varies  markedly  because 
of  changes  in  its  thickness.  The  actual  temporature  diffcrcnc.:  tlirough 
the  thcrmocline  is  a  relatively  stable  feature  and  is  uh-'rafore  predic¬ 
table  jll]  . 
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loy  n,  ;.~rl  'i  ■,..  ..  V.-?  ^  ir.Vctjc  co\’J-c’  fcr).  :iyd  ,;roya3'',to  di’rin,';  t!:.:;  nui.sior 


at  the  level  of  tlij  eeaconal  t.'ier:.iocllno  at  Ocean  Station  P^pa.  Vj'ith 
one-hoiirly  obc. ervatione  available  for  analysis,  the  aiinimn  neried  that 
conld  be  Invortipated  is  taro  hoairs. 


The  'Pukey  spcctrim  analysis,  proyraairied  for 
i.'as  used  to  analyse  each  tine  series.  The  min 
eredoriinant  periods  and  •iscoclated  anplitudes  pr 


the  3 DC  lh04  ccr.iputcr, 
object  y'as  to  find  any 
esont  in  tl:c  IXD  oscil¬ 


lations,  The  Tnkey  analysis  cctimtos  the  enor^y-densitp'  as 
of  i^criod  in  a  yiven  series.  This  is  done  by  conpietiny  tlio 
tion  function  for  a  nur.ber  of  yiven  lags  in  t!'.o  series;  then 


a  function 
autc-corrcla- 


a  snoothed 


spoctro,!  ostir.ato  is  obtain'.jd 


fron  the  Fourier  ccsinc  transform,  of  the 


auto-correlation  function. 

Initially,  a  corroloyran  ’.,'as  plotted  for  both  time  series  (figure  11), 
The  corrclocra::i  is  a  plot  of  the  non-nornaliacd  auto-corrclation  coeffi¬ 
cients  as  a  function  of  the  lag  (ic)  ,  Inspection  of  this  corrGlogr<ar.i  appears 
to  indicate  tliat  tlic  spectrum  of  series  1  has  a  relatively  Large  energy- 
density  inxinuji  near  the  12-hour  period  and  a  secondary  maxii.juri  vrith  a 
peak  near  six  hoars;  and  that  the  spcctnnn  of  series  2  has  only  one  naxinun 
near  12  hours. 


Koxt,  the  energy-density  spectrum  of  both  tine  scries  v?as  obtained 
from  the  Tulcey  analysis.  Figure  12  i  s  a  plot  of  the  energy-density  spectrum 
for  series  1.  A  maxinum  lag  (n)  of  24  hoixrs  \;as  used,  yhich  alloyed  the 
spectrum  to  include  a  range  of  periods  from  tvo  hours  through  4S  hours. 

The  spcctryn  helps  confirm,  the  indications  of  the  corrclogram.  A  rrucimura 
of  energy  of  the  variance  of  the  liT.D  is  shoxm  to  be  produced  by  fluctua¬ 
tions  with  periods  on  the  order  of  12  hours.  Also,  a  large  energy-density 
is  indicated  for  periods  around  six  hours.  There  is  an  energy  gap  for  the 
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periods  eround  eight  hours,  and  the  energy  for  periods  less  than  five  hours 
drops  off  rapidly  to  oji  insignificant  ai.iount. 

The  total  energy  is  a  known  quantity,  equal  to  twice  the  variance. 

By  measuring  total  area  under  the  spectrum  curve  and  the  area  under  each 
frequency  peak  bctvjcon  half-energy  points,  equivalent  amplitudes  for  v;aves 

a 

\7ith  periods  near  12  hours  and  six  hours  vroro  computed.  These  computed 
amplitudes  arc  18  feet  and  12  feet  for  the  12-hour  and  the  6-hour  periods, 
respectively. 

The  energy-density  spectrim  for  scries  2  is  shovm  in  figure  13«  A 
narciraun  lag  (n)  of  30  hours  \-7as  used,  which  allovied  inspection  of  periods 
ranging  from  tv7o  hours  through  60  hours,  A  relatively  large  energy  maxi¬ 
mum  is  again  indicated  for  periods  around  12  hours.  The  curve  then  falls 
off  rapidly,  but  shov;s  a  slight  maximum  near  the  6-hour  range  of  periods. 

A  single  sine  V7avc  with  energy  equivalent  to  that  under  the  12-hour  pealc 
would  have  an  amplitude  of  nine  feet. 

This  analysis  reveals  three  features  of  the  thermal  stiaicttu'c  at 
Ocean  Station  Pana  during  the  tine  intcrva.1  investigated, 

(l)  Vertical  oscillations  of  the  MLD  are  present  at  the  thcrmoclinc 
during  the  cooling  season.  Tne  magnitudes  of  these  oscillations  are 
similar  in  all  data,  \;hcther  they  be  from  hourly,  txiicc-daily,  or  daily 
obsor'/ations. 

^2)  Although  the  time  series  available  \7as  much  too  short  on  which 
to  base  strong  conclusions,  the  analysis  did  reveal  that  the  raajor  energy 
peaks  occurred  \Jith  a  period  of  approximately  12  hours  vrith  an  average 
amplitude  of  13.5  feet. 

(3)  The  value  of  the  standard  error  of  estimate  resulting  from  the 

BUD  07  regression  equations  is  very  close  to  the  height  of  the  postulated 
12-hour  wave  (that  is,  27  feet).  Therefore,  it  is  concluded  that  there 
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ic  a  direct  relation  botuoon  this  error  and  the  presence  of  internal 


\:avcs  during;  the  period  of  investigation. 
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7,  Conclunionc  and  Rocoiicacndationn 

As  a  result  of  this  study  of  the  influences  on  the  TiLD  during  the 
cooling  season,  the  follo\;ing  conclusions  can  be  reachea. 

(1)  Thernohalinc  convection  is  the  important  process  for  the  deepen¬ 
ing  of  the  ML'D  in  the  autumn  season. 

(2)  Dei  ping  the  period  of  study,  the  M.D  is  too  deep  to  bo  affected 
by  v;avc  stirring  or  diurnal  trends. 

(3)  Th2  tine  lag  bct\;ccn  surface  cooling  and  deepening  of  the  ILD 

by  the  convective  process  is  snail;  this  caiggests  that  once  the  cooling 
season  begins,  convective  circulation  is  essentially  continuous.  This 
circulation  nay  be  maintained  during  periods  of  near-zero  heat  across 

the  air-sea  interface  by  wind  energy,  as  suggested  by  I.aFond  [«]  ;  hov:- 
cver,  \mtil  an  injection  of  potential  energy  in  the  form  of  rapid  surface 
cooling,  the  M.D  uill  remain  relatively  static.  If  the  increased  surface 
density  caused  by  sui'face  cooling  is  sufficient,  the  convective  circula¬ 
tion  intensifies  and  drives  the  MLD  deeper. 

(4)  A  major  portion  of  the  MjD  variability  appears  to  be  associated 
■v.'ith  internal  \.'aves.  The  analysis  of  internal  v/aves  revealed  that  the 
raajor  energy  pealcs  occurred  with  a  period  of  approcinatcly  12  hours  v;ith 
an  average  amplitude  of  13.5  feet, 

Karly  in  1953,  at  a  Navy- sponsored  thornoclinc  conference  attended 
by  20  of  this  country’s  most  prominent  oceanographers  and  meteorologists, 
the  follo\.’ing  question  uas  asked  his  colleagues  by  Dr.  T.  F.  Malone;  "How 
accurate  a  prediction  of  the  MLD  and  its  hourly  fluctxiations  is  required 
by  the  military?"  [19]  This  remains  a  vital  question, 

Tf  the  military  is  \/illing  to  accept  forecasts  of  the  MLD  within 
ranges  of  20  or  30  feet,  then  existing  methods  of  prediction  will  suffice. 
Any  increase  in  accuracy  over  existing  schemes  for  predicting  the  MLD  at 
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a  n'’,i’ticul’r  lioiu'  on  a  narticular  day  can  bo  acconpliohcd  in  aroac  vrhorc 
internal  viavcc  aro  inportant  only  after  prediction  of  tho  internal  \:avc 
pattern  is  nossible.  This  is  a  conplicatcd  problen. 

If  a  sufficiently  lon^,  detailed  tir.ic  series  of  BT's  could  bo  mdc 
available,  intornal  v;avcs  could  be  filtered  out  by  use  of  a  technique  such 
as  proposed  by  Linnettc  [u].  This  uould  leave  a  tine  series  of  the  M.D 
that  is  due  entirely  to  convective  miring  or  other  causes,  thus  enabling 
a  mch  better  quantitative  analysis  of  tho  convective  process.  This  is 
suggested  as  the  next  logical  extension  of  the  research  reported  in  this 
paper. 
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Figure  PoGition-Indicating  Grid  for  Ocean  Weather 
Station  Papa,  uith  a  Ilercator  Projection  of 
a  Lao-tude  and  Longitude  Grid  Superinposod 
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Fi^xiro  5»“'*  CroDK—Spec  -  XD  and  Sova-Svirfaco - 

ToiiporatuTG  Tine  Series  (llovoraber^  19^S) 
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Figure  7.  Average  Daily  MLD 

(October  through  December,  1956) 
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Figure  11.  Correlogram  for  One-Hourly  Time  Series 
(October,  1957,  and  Novembfer ,  1958) 
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Table  I 
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T:'J)1e  II 


02005 
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.039 

(P)^5 

-.060 

-.004 

-.101 

.041 

GPM5 

.118 

-.463 
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.091 
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020C5  and  170CS 

1957 
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1956 
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1958 
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-.489 
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-.020 

.017 
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.092 
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.290 

.202 

.1/a 

.263 

.238 
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.327 

.148 
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.276 

.004 

.231 

.307 
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.047 

.360 
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.487 

-.058 

-.127 

.581 

-.060 

-.649 

-.843 

-.899 

-.726 

-.819 

1 

1 

-.314 

-.505 

-.502 

-.277 

.383 

.352 
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-.143 
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-.115 

.525 

-.389 

-.192 

.583 

-.245 

.170 

-.342 

-.107 

-.027 

-.090 

.327 

-.215 

.047 
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-.419 

.143 
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-.050 
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.028 
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Table  III 


Soa-S’oi’f  ac  o-  reMpr-ratui’G 
Lar  in  Koiirs 

0 

2 

5 

0 

o 

11 

1/t 

17 

20 

23 


Coi-rolatioa  Coofficicn 
_ witli  the  IliD 

-.04975 

-.88409 

-.87828 

-.37576 

-.37173 

-.86479 

-.86829 

-.86125 

-.85189 


Sinplo  Linoar-Corrclation  Coefficients 
for  Sor^S’arf'LCo  Te.  loeratures  and  the  IXD 
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Table  IV 


Linear-correlation  coofficiont 
(for  To_2)  -.S84 

Standard  cri*or  of  the  estimte  30.33 

Range  of  residuals  121.2 

Regression  equation;  =  1152.94  -  20.24  (Tj,_2) 


Regression  Results  for  One  Independent  Variable 
(17003,  1956) 
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Table  V 


1956 

0200a 

1957 

1958 

1956 

1700S 

1957 

1958 

02002  and  17002 

1956  1957  1958 

(Ty  and 

T^,;_2  included,  T 

a-Ty  and  (Tg^-Ty). 

.2  onitted) 

.950 

.941 

.864 

.915 

.831 

.915 

.909 

.874 

.863 

26.5 

14.4 

26.7 

31.4 

18.0 

19.4 

29.9 

15.8 

21.0 

80.4 

45.4 

29.9 

104.8 

49.6 

59.8 

118.8 

57.1 

102.5 

210.1 

179.4 

273.3 

212.7 

170.5 

279.6 

211.6 

174.8 

276.8 

67.3 

32.4 

37.0 

64.9 

25.0 

36.27 

65.7 

29.0 

36.5 

and  (T^- 

•lvr)_2  included, 

Ty  and  Ty. 

__2  onitted) 

.950 

.917 

.910 

26.5 

30.9 

29.8 

80.4 

104.4 

121.8 

210.1 

212.7 

211.6 

67.3 

64.9 

65.7 

Ihltiple-Regrossion  Recults  for  18  Independent  Variables 


Table  VI 


I'iultiplc-sorrelation  coefficient 


(for  18  independent  variables)  .896 

Linear-correlation  coeffici ''nt 

(for  Tg  and  Tg_2)  -.814 

Standard  error  of  the  estimte  22.15 

IlanfTe  of  residuals  72.0 

Standard  deviation  43. B 


Thi'oe-Year  Avoi'af^es  for  Regression  Results 
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.4PEi:ni)ix 


Heat  Transnort  in  t’no  Sea  OJ^d  ?.xchan)"e  acronr.  Air-Soa  Interface 

A.  Thernohaline  convection.  The  vertical  displacencnt  of  orall  oar- 
cels  of  nater  occurs  U’hon  a  srall  part  of  n  v/ater  nass  is  heavier  than  the 
i,;ator  underneath  it.  To  restore  this  unstable  condition  to  eqTiilibrivm, 
the  heavier  water  tends  to  sink  \/hile  the  lighter  v/ater  rises.  Associated 
\/ith  these  forced  vertical  novements  of  snail  parcels  there  is  also  a  trans¬ 
port  of  the  chara.cteriEtic  properties  of  oca  water  in  a  vortical  direction 
which  leads  to  an  equalisation  of  any  vertical  differences  in  these  propci*- 
tlco  'Jhich  my  be  present.  It  is  believed  that  convective  ni:cing  is  as 
iraportant  an  agent  of  heat  transfer  in  the  ocean  as  wind  stirring  during 
the  cooling  season  ^2,  12.  2(^  , 

An  initial  increase  in  the  density  of  smll  particles  at  the  corface 
accorapanics  an  increase  in  salinity  (due  to  evaporation,  or  to  the  form- 
tion  of  ice)  or  a  decrease  in  temperature.  One  or  more  of  these  in  combi¬ 
nation  my  be  involved  in  thernohaline  convection.  In  lower  latitudes, 
v/here  there  are  only  smll  variations  in  the  temperature,  the  heat  loss 
is  outweighed  by  the  effect  of  evaporation;  in  polar  regions,  in  addition 
to  radiation  and  evaporation  the  increase  in  salinity  due  to  formtion  of 
ice  is  also  effective.  In  temperate  latitudes  the  heat  loss  by  radiation 
is  the  decisive  factor. 

Convective  sensible  heat  exchange  (conduction)  between  air  and  sea 
surface  has  been  studied  by  Kuhbrodt  [A]  for  the  North  Atlantic.  He  has 
shoim  the  convectional  heat  flux  to  be  anproximatcly  20  gm  cal/cn^  per  day 
during  the  period  of  the  year  when  the  air  is  cooler  than  the  sea  surface. 
This  value,  oven  as  a  rough  estimte,  is  too  large  to  be  neglected  in  the 
heat  budget  of  the  ocean.  If  the  surface  of  the  sea  is  warmer  than  the 
atmosphere,  the  air  is  heated  at  the  interface  and  the  vcrt3c.-xl  stratifi- 
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cation  of  the  air  becomes  unstable;  the  air  bccor.GS  tui'bulcnt  -./ith  a  con¬ 
stant  replacement  by  cool  air  at  tliC  sea  surface  and  the  vertical  heat 
transport  becomes  large. 

Values  given  by  Defant  for  the  heat  loss  in  \inito  of  fjr  jal/cm^ 
per  day  for  50N  latitude  are  116  units  duo  to  effective  back  radiation, 

78  units  duo  to  evaporation,  and  20  ixnits  due  to  convection  of  sensible 
heat  (conduction)  as  a  yearly  average.  Certainly,  during  the  cooling 
season  heat  loss  due  to  conduction  uould  be  even  more  significant.  These 
values  are  derived,  assuming  that  the  heat  exchange  through  the  ocean  sur¬ 
face  occurs  independently  for  each  separate  latitude  belt.  Therefore,  no 
meridional  heat  exchange  (by  ocean  currents  and  by  horizontal  nixing)  \;ao 
considered, 

B.  Dynamic  convection  (forced  vertical  nixing  due  to  wind  and  wave 
action) .  This  process  has  been  v:ell  studied  and  considered  to  be  the 
dominant  process  of  vertical  nixing  in  the  heating  season.  Laevastu  (l?^ 
has  proposed  an  empirical  formula  for  determining  the  MLD  by  significant 
v;ave  height  alone. 

C.  Dynamic  convection  due  to  forced  vertical  nixing  of  ocean  currents 
can  result  in  convective  nixing,  and  influences  the  ME,  Tliis 
process  \/as  not  considered  because  of  the  difficulty  of  measurement;  hov/- 
evor,  since  the  Station  is  located  in  the  center  of  the  Alaskan  Gyral, 
ocean— current  effects  are  probably  smll, 

D.  Thermal  conductivity  can  take  place  ’..'hen  a  ^'ertical  temperature 
gradient  exists  in  the  ocea:' ,  Heat  is  transferred  by  molecular  heat  con¬ 
duction  processes.  Defant  |4j  concludes  that  this  process  is  insignificant 
in  oce.anogranhic  investigations  due  to  the  extremely  long  period  involved 
in  the  conduction  process.  Therefore,  this  process  v;ill  not  be  considered 
in  this  study. 
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^ntrai!i:nont  d~'o  to  turbulencs  at  tl’o  Ycnoclinc  ♦  Durin;;;  the 
coolin;'  acason  in  the  open  ocean,  a  v/ell-clefined  ioothornal  surface  layer 
exists  v.’hich  e::tcnds  dovm.'ard  to  an  interface  of  rapidly  decreasinp;  ten- 
)3oratui'c,  Investi';;ationG  of  the  subarctic  water  indicate  an  isohaline 
condition  to  a  depth  of  331  feet,  so  the  effect  of  salinity  on  density 
can  be  ignored  [23J ,  and  thus  mo  my  consider  the  therriodine  as  a  pycno- 
cline.  Because  less  dense  wari.ier  v/ater  overlies  dense  cold  \:ater,  v/e  have 
a  stable  condition  mch  like  a  surface  inversion  in  the  atr.iospliere,  When 
surface  cooling  takes  place,  thereby  crea.ting  an  unstable  condition,  free 
convection  can  occui’.  Cronwcll|”3j  has  denonstrated  that  this  free  convec¬ 
tion  in  a  nodel  causes  a  turbulent  exchange  across  the  thernocline  vfhich 
is  ono-way.  The  fluid  particles  ;/hich  wove  upvfard  from  the  region  beloi; 
the  thernocline  are  rapidly  deforned  and  laixed  tiiroughout  the  upper  layers; 
fluid  particles  ;.'hich  nove  do;:mjai'd  tovfard  the  quiet  layer  below  the  ILD 
are  buoyed  upward  intact.  Thus,  the  upper  layer  increases  in  thickness, 
but  decreases  in  temperature, it  the  expense  of  the  lo\;er  layer.  This  sane 
phenonenon  is  observed  when  the  turbulence  is  induced  by  v/ave  action  or  in¬ 
ternal  wave  action  [l]  . 
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